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THE PULSAR B2224+65 AND ITS JETS: A TWO EPOCH X-RAY ANALYSIS 
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ABSTRACT 

We present an X-ray morphological and spectroscopic study of the pulsar B2224+65 and its apparent 
jet-like X-ray features based on two epoch Chandra observations. The main X-ray feature, which 
shows a large directional offset from the ram-pressure confined pulsar wind nebula (Guitar Nebula), 
is broader in apparent width and shows evidence for spectral hardening (at 95 percent confidence) in 
the second epoch compared to the first. Furthermore, the sharp leading edge of the feature is found to 
have a proper motion consistent with that of the pulsar (~180 mas yr~^). The combined data set also 
provides evidence for the presence of a counter feature, albeit substantially fainter and shorter than the 
main one. Additional spectral trends along the major and minor axes of the feature are only marginally 
detected in the two epoch data, including softening counter to the direction of proper motion. Possible 
explanations for the X-ray features include diffuse energetic particles being confined by an organized 
ambient magnetic field as well as a simple ballistic jet interpretation; however, the former may have 
difficulty in explaining observed spectral trends between epochs and along the feature's major axis 
whereas the latter may struggle to elucidate its linearity. Given the low counting statistics available in 
the two epoch observations, it remains difficult to determine a physical production scenario for these 
enigmatic X-ray emitting features with any certainty. 
Subject headings: 



1. INTRODUCTION 

As pulsars age, they lose their spin-down energy 
through both radiation and relativistic particle ejection. 
One may approximate the particle ejection in two forms: 
anisotropic winds from the pulsar surface and bi-polar 
jets. Previous X-ray studies have shown the wind com- 
ponent to be primarily equatorial and axisymmetric, pro- 
ducing observed tori in pulsar wind nebulae (PWNe), 
whereas the jet outflows are only observed in a few, 
though well studied, young pulsars (e.g. the Crab an d 
Vela pulsars, IWeisskopf et"alll2000t iPavlov et al]|2003[ ). 
Only a few nearby and relatively old pulsars have been 
shown to have their PWN e, mostly enhanced by ram- 
pressure confinement (e.g. , IWang. Li fc Begelmanlll993t 
iKargaltsev fc Pavlov! I2008L and references therein). For 
a PWN or similar pulsar ejection feature to be observed 
in X-rays, the synchrotron-emitting particles must have 
energies on order of 10^ TeV, achievable in shocks of 
the pulsar ejecta as it encounters the local interstellar 
medium (ISM). X-ray investigations of pulsar outflows 
are therefore crucial in understanding their origin and 
particle acceleration, which can in turn provide insight 
into the origin of other, more exotic and distant, outflows 
(e.g. from Active Galactic Nuclei). 

A particularly striking case is that of the pulsar 
B2224-f65 and its X-ray-emitting "jet". B2224+65 
is one of the fastest pulsars known with a ra- 
dio proper motion of M = 182 ± 3 mas yr~^ 
([Harrison. Lvne fc AndersonI Il993l ): corresponding to a 
transverse velocity of 864 Dkpc km s~^, where Dkpc is 
the distance to B22244-65 in uni ts of kpc and is un- 
certain in the range of ~ 1-2 ()Chatteriee fc CordesI 
|2004[ ). The rapid motion of the pulsar produces a bow 
shock neb ula visible in H a (Cordes. Rom ani fc LundgrenI 
[1993 : Chatteriee fc Cordes 2002 . 2004 ) which has been 
dubbed the "Guitar Nebula" due to its peculiar mor- 



phology. Aside from its high proper motion, B2224-I-65 
behaves like a normal radio pulsar, with a modest en- 
ergy spin-down rate oi E = lO^'^'^ ergs s~^ and period 
P = 0.68 s. X-ray studies of B2224+65 do not reveal sig- 
nificant emission concurrent wi th the Ha (however, see 
iRomani. Cordes fc Yadigaroglul[T997[ ) but rather an un- 
usual extended linear feature apparently stemming from 
B2224H-65 and offset from the pulsa r 's direction of mo- 
tion bv ^118° jHui fc Becked ^OOTI: iWong et all [20031 : 
IZavlin fc Pavlovll2004l ). The feature, extending ~2' away 
from the pulsar, was found to have an apparent non- 
thermal spectrum though the details as to its origin have 
thus far remained a mystery. Such a distinct, linear X- 
ray emitting feature has not been identified anywhere else 
in the Galaxy with the possible excep tion of the Galac- 
tic Center (e.g. I Johnson. Dong, fc VV^ ng 2009, and refer- 
ences therein). It was suggested bv [Kargaltsev fc Pavlov! 
(j2008) that the extended feature may not be associated 
with B2224-I-65 but rather with a nearby point source, 
although no physical production scenario was developed 
for the feature. Throug h theoret i cal m odeling of the lin- 
ear X-ray emission by .Bandiera! ([2008! ) , it was proposed 
that the feature could be produced through highly ener- 
getic particles that are accelerated at, and subsequently 
escape, the pulsar wind terminal shock but are ultimately 
confined by a large-scale, organized interstellar magnetic 
field. However, th ere are a numbe r of observational re- 
quirements to the !Bandiera[ ([2008f ) model that have yet 
to be fulfilled, primarily a proper motion of the feature 
consistent with B2224+65. The uncertain association 
with the pulsar, combined with low photon statistics, 
has made it difficult to establish the physical process re- 
sponsible for the feature. 

In this study, we examine two epoch data taken 
by the Chandra X-ray Observatory to determine any 
association between the extended, linear feature and 
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TABLE 1 
Chandra ACIS-S Observations 



Chandra 


R.A. (J2000) 


Dec. (J2000) 


Exposure 


Roll Angle 


OBS Date 


OBS. # 


(h m s) 


(degree arcmin arcsec) 


(s) 


(degree) 


(yyyy-mm-dd) 


755 


22 25 48.439 


+65 35 08.09 


47865 


236.8 


2000-10-21 


6691 


22 25 51.852 


+65 35 39.19 


9968 


182.7 


2006-08-29 


7400 


22 25 51.428 


+65 35 35.13 


35839 


219.3 


2006-10-06 



Note. — The exposure represents the live time (dead time corrected) of cleaned data. 
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Fig. 1. — Merged epoch Chandra intensity map of B2224+65 in the 0.3-7.0 keV energy band. The image has been smoothed using 
an adaptive Gaussian filter (FWHM<1") and is presented with logarithmic scaling. Overlaid on the intensity map is the direction of 
B2224+65's proper motion along with detected X-ray sources. Also included are the regions used for spectral extraction of the pulsar and 
jet (solid line), the local background region (the area within the dotted line excluding sources and the dotted rectangle), and the regions 
used in constructing surface brightness profiles along the major (dot-dash line) and minor (dashed line) axis. These regions all have position 
angles of '^24.7° as determined from the bright edge of the jet within '^50" from B2224+65. 

B2224+65 as well as to place constraints on possible X- 
ray-emitting particle production scenarios. Throughout 
this manuscript, we will refer to the feature as the jet or 
jet-like feature based purely on its apparent morphology, 
implying no physical interpretation given its still uncer- 
tain origin. From the two epoch data, we first correct for 
any astromctric offset between observations and verify 
the motion of B2224+65. We then briefly examine the 
two epoch and combined spectroscopic information for 
the pulsar to see if there may have been any evolution 
in its spectrum. A morphological analysis of the jet will 
serve to detect any shift between epochs, providing evi- 
dence for (or against) an association between the jet and 
pulsar. Spectral analysis of the jet between epochs and 



across the major/minor axes may aid in identifying pos- 
sible spectral trends, which can assist in determining the 
physical models behind the X-ray emission. F inally, we 
briefly discuss our findings in the context of the lBandieral 
(|2008( ) model and a simple ballistic jet model similar to 
that used to explain AGN jet emission. 

2. OBSERVATIONS AND DATA REDUCTION 

B2224+65 was observed by the Chandra X-ray Ob- 
servatory over two separate epochs using the Advanced 
CCD Imaging Spectrometer spectroscopic array (ACIS- 
S). The first epoch consists of Obs. ID 755 while the sec- 
ond epoch is composed of two observations, Obs. ID 6691 
and 7400 (the relevant observational parameters for each 



Fig. 2. — Intensity maps of B2224+65 and the jet with the source region from Fig. 1 overlaid; 0.3-7.0 keV images in the first (a) and 
second (b) epochs as well as the merged images in the soft (0.3-1.5 keV; c) and hard (1.5-7.0 keV; d) bands. All images are adaptively 
smoothed similar to Fig. 1 and are shown here using the same relative logarithmic scale with the minimum set by the median image value 
and the maximum set a factor of 10^ times higher. 



Obs. ID are given in Table 1). Each of the individual ob- 
servations was first reprocessed using the standard Chan- 
dra Interactive Analysis of Observations (ciAO) routines 
(version 3.4.0). After reprocessing, Obs. ID 6691 and 
7400 were combined to produce a second epoch event file. 
A merged epoch event file and X-ray image of B2224+65 
(Figs. 1,2) was then produced by combining Obs. ID 755 
and the second epoch event file after astrometric cor- 
rections have been applied (see below). Source lists for 
each of the two e poch event fil es were constructed us- 
ing the method of IWansj (|2004f ). with a false detection 
probability of 10~^, to be used in correcting the astrome- 



try between epochs and producing source-removed event 
maps for background estimation. 

Before we can properly analyze the two epoch data 
sets, we must correct for any possible astrometric offset 
that may have been produced through pointing differ- 
ences between epochs. From our X-ray detected source 
lists, we extract the positions of 19 sources common to 
both epochs with count rates >0.2 counts ks~^ and off- 
axis distances <4'. Limiting the off-axis distance of the 
sources serves to minimize the effects of positional errors 
due to PSF variation across the CCD. The RA and DEC 
offsets, as well as a rotation angle in image alignment, 
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are then derived from the RA and DEC differences of 
the common sources, taking into account their positional 
uncertainty determined from maximum hkehhood cen- 
troiding (with la errors of ^0.1-0.9"). From our initial 
calculation, we find that the rotation between images is 
negligible (0.02±0.05°, Icr uncertainty) and subsequently 
recompute the corrections assuming no rotation, result- 
ing in ARA and ADEC corrections of 0.24±0.09" and 
0.18±0.09", respectively. Despite exclusion of the rota- 
tion in image alignment, it is unlikely to influence our 
results, as we limit our analysis to only the area near the 
telescope aim-point. These astrometric corrections have 
been applied only to the astrometry in the fits header 
of the second epoch event file for merged epoch image 
production (Figs. 1,2) and subsequent analyses. 

Due to the diffuse nature of the jet, we adopt the 
method of local background estimation to determine the 
background spectra. To obtain a representative back- 
ground for the pulsar and jet regions, we extract the 
spectrum from a ~2' circular region, given by the dotted 
circle in Fig. 1, while excluding detected point sources 
and the jet. Since there may be additional, unresolved 
emission from the jet extending beyond its source region, 
we exclude the area within the dotted rectangle of Fig. 1 
when extracting the background. This same area is also 
used for determining the sky background level when con- 
structing surface brightness profiles along the major and 
minor axes of the jet (§ 3.2). 

3. ANALYSIS AND RESULTS 

3.L Pulsar B2224+65: Proper Motion and Spectrum 

Due to the time gap between the observations of 
B2224-f65, we must first verify the pulsar's proper mo- 
tion {fiRA = 144±3 and (Idec = 112 ± 3 or ^ = 182 ± 3 
mas yr^^ with 90 perc ent confidence inte rvals from the 
radio sample of Harrison. Lvne fc Anders on f993). Ex- 
tracting B2224-|-65's position from our X-ray source lists 
in both epochs and applying the astrometric correc- 
tions, we find that the pulsar's position has changed by 
0.91±0.24 arcsec in RA and 0.44±0.24 arcsec in DEC, 
taking into consideration the relative astrometric uncer- 
tainties, or by 1.01±0.34 arcsec over the '^6 year gap 
between epochs. This relates to a proper motion of 
fijiA ~ f53±40 mas yr~^ and /idec ~ 74 ±40 mas yr~^ 
or absolute proper motion of « 170 ± 57 mas yr"-'^ (Icr 
errors), in agreement with the radio proper motion. 

While we do not expect variation in the spectrum of 
B2224-I-65 itself, a significant portion of the observed 
emission may arise from the shocked wind material, 
which is dependent on the pulsar wind and local ISM 
properties, and could therefore change between epochs. 
Any variation, or lack thereof, in the emission can then 
be compared with the spectral modeling of the jet. To 
ensure that the pulsar's proper motion will not influ- 
ence the spectral modeling, we extract the spectra from 
the individual epoch event files using a small circular re- 
gion (radius of 3" at RA=22h25m53s, DEC=-H65°35'35" 
(J2000), see Fig. 1). The spectra are then adaptively 
binned such that each bin has a signal-to- noise (S/N) ra- 
tio of >2.5 where the signal is the net counts above the 
local background and the noise contains all counts. This 
provides 8(9) bins for the first(second) epoch spectrum 
with ^10 counts per bin due to the small contribution 
from the local background over the pulsar region. 



To model the two epoch pulsar spectra, we apply a 
joint fit, absorbed power law in XSPEC (version 11.3.1) 
where the foreground X-ray-absorbing column density, 
Nh, should not change significantly with time and is 
thus set as a common parameter between epochs. The 
epoch averaged properties of the pulsar are derived by 
performing a similar joint fit in which all model param- 
eters are set to be the same in both epochs (i.e. power 
law index and normalization). Alternatively, one may 
also combine the two epoch spectra using the ADDSPEC 
routine in ftools and apply the absorbed power law 
model, or directly extract the spectrum from the merged 
event file. However, the latter method does not preserve 
calibration information from each epoch whereas the for- 
mer may give misleading results if the detector response 
changes between spectra and/or if there is any spatial 
variation in the spectra. Nevertheless, we find that each 
method produces consistent results with the simple two 
epoch joint fit. We report here only the spectral param- 
eters as given by the joint fit modeling as it is likely the 
least sensitive to any systematic calibration uncertain- 
ties. Comparison of our spectral fitting results (with 90 
percent confidence contours given in Tabl e 2 and Fig. 3) 
of B2224-I-65 with the previous work of iHui fc Becked 
(j2007D (Icr errors based on first epoch data alone) sug- 
gests that the observed emission from the pulsar has not 
had any significant evolution between epochs. 

3.2. Jet-Like Feature: Astrometry and Spectrum 

From the smoothed X-ray intensity map of B2224+65 
(Fig. 1), the jet appears disconnected from both 
B2224-f 65 and the nearby point source (labeled here as 
CXO J222552.0-h653548) and thus does not clearly sup- 
port a direct asso ciation with either source (c ontrary to 
the suggestion of iKargaltsev fc Pavlov! [20081 ) . The jet 
remains fairly linear with a sharp edge along the direc- 
tion of B2224-|-65's proper motion, extending a length of 
'^2.4'away from the pulsar. Other morphological proper- 
ties include apparent "fanning" with increasing distance 
from the pulsar and a minor region of X-ray enhance- 
ment on the opposite side of the pulsar. Furthermore, the 
two epoch images (Fig. 2a, b) indicate an increase in the 
width of the jet across epochs with little to no bending in 
the leading edge (<1° within -2' away from B2224-h65) 
while the soft and hard band images (Fig. 2c, d) show evi- 
dence for harder X-ray emission farther along the (epoch 
averaged) jet. In the following analysis, we assume the 
case of no bending for simplicity. The presence of any 
bending will serve primarily to dilute the significance of 
any determined edge location along the minor axis and 
will be propagated through our analysis as a systematic 
uncertainty. 

For our morphological analysis of the jet, particu- 
larly for studying its proper motion and the signifi- 
cance of the extended emission regions, we construct 
surface brightness profiles along its major and mi- 
nor axes. These profiles result from 100x25 arcsec^ 
and 85x50 arcsec^ rectangular regions centered at 
RA=(22h25m52s, 22h25m45s) and DEC=(+65°35'35", 
-|-65°35'57") (J2000) for the major and minor axis, re- 
spectively (given as the dot-dash and dashed line regions 
in Fig. 1). These regions were selected by hand to max- 
imize the signal for the relevant profile using the epoch 
averaged pulsar position as the primary reference point 
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TABLE 2 
Spectral Fitting of B2224+65 



Epoch 


Count Rate (0.3-7.0 keV) 
10-3 cts/s 


10^2 cm-2 


r 


Flux (0.3-7.0 keV) 
lO"^* ergs cm~^s~^ 


xVd.o.f. 


1 
2 

Merged 


1.73±0.19 
2.00±0.21 
1.86±0.14 


<0.09 
<0.09 


l.DD_Q 26 


„ +0.44 
^•'^-0.64 
1 gi+0.50 

, OO+0.22 
t-«3_o 3g 


8.94/12 
9.47/14 



Note. — xspec spectral fitting parameters for B2224-I-65. The column densities listed are the 
90 percent confidence upper limits given by the absorbed power law fits. Errors are given to 90 
percent confidence. 
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Fig. 3. — 68 percent (dashed line) and 90 percent (solid line ) 
confidence contours for the epoch averaged pulsar spectral fit for 
Njj versus F and for F versus normalization. The crosses in the 
plots mark the best-fit parameter values from the joint fit. 



(§ 3.1). For the calculation of the profiles, we bin the 
profiles with a S/N ratio of >4; here the signal is the net 
counts above the instrument background, as measured 
from the ACIS stowed data, while the noise includes all 
counts. The use of the instrument background, instead 
of the total sky background, enables the application of 
the adaptive binning even to regions with negligible con- 
tribution from the jet. For each bin, the X-ray inten- 
sity is then simply defined as (counts-instrument back- 
ground)/exposure. After binning, the sky background 
levels, determined from the local background region, 
were subtracted to produce the surface brightness pro- 
files shown in Fig. 4. The adopted binning was chosen 
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Fig. 4. — X-ray intensity above the sky level along the major (a) 
and minor (b) axes, adaptively binned to a S/N of >4 per bin. For 
reference, the sky background levels for each profile are 51.9±0.9, 
74.5±1.2, and 62.6±0.7 xlO"* counts s"! arcmin'^ for the first 
epoch minor axis, second epoch minor axis and combined epoch 
major axis profiles, respectively. Distances are calculated relative 
to the epoch averaged pulsar position. Errors are presented at the 
1(T confidence level. 

simply for ease of viewing; the exact choice does not in- 
fluence our analysis of the jet edge or extended emission 
regions (see below). 

Examining Fig. 4a, there is evidence for a weak counter 
feature, opposite in direction to the primary jet-like fea- 
ture. This counter feature, extending '^20" away from 
B2224-f 65, is considerably fainter than the main jet-like 
feature though it still appears significant compared to 
the local X-ray background. Due to the faintness of the 
counter feature, the counting statistics are too low in in- 
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Fig. 5. — Histogram distributions of edge locations in the first 
and second epochs. The distributions result from bootstrap re- 
sampling of count profiles along the minor axis, taken from a fine- 
binned count map, and convolving the realized profiles with the 
chosen filter. The x-axis gives edge locations relative to the epoch 
averaged pulsar position in arcsec. 

dividual epoch observations for a positive detection. We 
therefore use the merged epoch event file, along with a 
small 20 x 20 arcsec^ box centered on the feature in ques- 
tion, to derive a background subtracted count rate of 
0.34 ± 0.10 counts ks~^. Comparatively, the local back- 
ground corresponding to the same region has a count rate 
of 0.68±0.09 counts ks~^, placing the counter feature at 
4 times the background uncertainty or at a detection 
significance of ^Aa. It is therefore highly unlikely that 
this feature results from Poisson variations in the X-ray 
background. 

In Fig. 4b, there appears to be a slight shift in the 
sharp left edge as well as an extension far to the right 
when comparing the minor axis surface brightness pro- 
files in the first and second epochs. Currently, we concern 
ourselves with the presence of the apparent shift which, 
if confirmed, may be compared to the observed proper 
motion of B2224-I-65, establishing a connection between 
the pulsar and jet-like feature. In order to determine the 
shift between epochs, we apply a simple one dimensional 
edge detection method to the sharp edges in the two 
epoch data. Our edge detection method involves con- 
structing finely binned intensity profiles taken from high 
resolution count maps (with ^0.12" pixels) and convolv- 
ing the profiles with a filter based on the work of ICannvl 

(fl986l) . which can be approximated hy xe ^ , with 
a filter scale of s=0.5", comparable to the ACIS resolu- 
tion. To determine the underlying distribution of edge 
positions in the two epoch data, and thus derive uncer- 
tainty estimates, we generate 5000 Poisson realizations 
of the high resolution count profiles and convolve each 
with the chosen filter, following the method of bootstrap 
resampling. By filter design, the edge position for each 
realization is determined by "centroiding" the peak in 
the surface brightness profile and therefore depends pri- 
marily on the counting statistics of the available data. 
The derived positions are thus insensitive to the exact 
choice of the filter scale or binning of the count profiles; 
we have verified this using a variety of profile binning 
and filter scales, all of which produce results consistent 
with those presented here. Based on the histogram dis- 



tributions of the resampling (Fig. 5), we find a shift be- 
tween the two epochs of 0.621;q'28 arcsec, with 90 percent 
confidence intervals. The reported errors on the shift 
are likely underestimates of the true uncertainty as they 
do not account for errors in our astrometric corrections, 
though they do include any possible systematic uncer- 
tainty from unresolved bending in the jet. Nevertheless, 
the shift, especially its direction, is consistent with the 
proper motion of B2224-I-65, correcting for its projection 
to the adopted minor axis of the jet (a factor of ^0.88), 
and thus supports an association between the pulsar and 
jet-like feature. 

As mentioned previously, there is evidence for excess 
extended emission along the minor axis of the jet in the 
second epoch compared to the first. This excess emission 
can also be seen in Fig. 2a, b through the apparent widen- 
ing of the jet-like feature over the two epochs. Examining 
the two epoch surface brightness profiles, there is little 
deviation between epochs, after correcting for the pre- 
viously determined shift, except for the excess emission 
region, located ~7-17" away from the left sharp edge. 
The presence of this region will not effect our previous 
edge determination due to its distance from the left sharp 
edge and our choice of filter. Following a method similar 
to the extension along the major axis, the region of ex- 
cess extended emission is found to have a count rate that 
is 0.28±0.11 counts s^^ higher in the second epoch than 
in the first, which has a total count rate (sky background 
included) of 0.64±0.06 counts s~^. Given the high de- 
tection significance of the excess emission in the second 
epoch (^5cr), it is very unlikely that it is artificial or sta- 
tistical in nature; although a physical explanation for the 
excess emission is perhaps even more enigmatic than the 
jet or counter feature. 

We now examine the spectrum of the jet both spatially 
and across epochs to determine if there may be evidence 
for any spectral evolution. In each epoch, we extract the 
spectrum of the jet using the same rectangular region 
(see Fig. 1,2) after removing B2224-I-65. The jet region 
is further segmented along its major and minor axes, la- 
beled according to their relative positions, such that the 
segments have approximately equal counting statistics. 
The individual spectra are again adaptively binned to a 
S/N ratio of >2.5, resulting in >25 counts per bin due to 
the higher background contribution in the jet region. In 
order to improve the constraints on model parameters in 
the spectral fits for the segments, we again apply joint fits 
in which the parameters are set to be the same in both 
epochs though allowed to vary over the individual seg- 
ments; only Nfj is set as a common parameter through- 
out, assuming that it is constant in time and space across 
the jet. We also include the epoch averaged pulsar spec- 
trum to place further constraints on Nh and the rele- 
vant Fs. The fitted parameters are given in Table 3 with 
the confidence contours for the epoch averaged jet spec- 
trum given in Fig. 6. Similar to Hui fc Be cker (2007), 
we find that the power law model provides a satisfac- 
tory fit (x^/d.o.f ~1.1) to the jet spectrum, though the 
possibility of alternative models can not be completely 
rejected based on test alone. It should be noted, 
however, that reasonable alternatives such as thermal 
Bremsstrahlung result in temperatur es on the order of 
>2 keV (see also lHui fc Beckeiir2007| ) and are likely un- 
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TABLE 3 

Spectral Fitting of the Jet-like Feature 



Epoch 


Segment 


CR (0.3-7.0 keV) 
10~^ cts/s 


10^2 cm-2 


F 


Flux (0.3-7.0 keV) 
10~^* ergs cm~^ s~^ 


X^/d.o.f. 










1 79+0.41 
^- '^-0.35 








Jet 


K 7q-|-n 4fi 


<0.12 


46.6/46 


2 


Jet 


5.d0±0.49 




1 ic:+0.43 
^■^■^-0.39 


■7 no + 1.99 

|-.3»_2 18 




Merged 


B2224:-\-bi> 


1.0D±U.14 




1 Qo+0.44 
^■'^^-0.37 


1 fiQ+0.30 










Epocii Averagc( 




Merged 


Jet 


5.70±0.34 


<0.15 


— 0.17 


5 58+°-'^'^ 


52.7/48 


Merged 


B2224-I-65 


1.86±0.14 




1 95+0.44 
i.yo_Q 21 


1 63+°'^« 










Major Axis 




Merged 


Head Segment 


3.61±0.25 


<0.13 


^■"^-0.34 


'-'•'-'^-0.84 


33.6/47 




Tail Segment 


1.96±0.21 




1 00+°-^^ 


q qc + O.SO 
1 13 






B2224-I-65 


1.86±0.14 




^■^'-'-0.37 


^•"0-0.40 










Minor Axis 




Merged 


Edge Segment 


2.64±0.22 


<0.06 


1 ns+0-43 


3 7fi+0.78 


31.9/45 




Taper Segment 


2.79±0.24 




1 so+°-^i 

J^-J'J_0.37 


q 7n+0.98 
•^■'"-0.76 
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— 68 percent (dashed line) and 90 percent (solid line) 
contours for the merged jet spectral fit similar to Fig. 3. 



realistic. From our spectral analysis, there is evidence for 
significant, to 95 percent confidence, spectral hardening 
in the jet between epochs. There is also indication for 
further spectral variation along the jet, including spectral 
hardening farther from the pulsar and softening counter 
to the pulsar proper motion, although these trends are 
only marginally detected (with ^ Icr significance). 

4. DISCUSSION 

Based on our two epoch X-ray analysis of B2224-I-65 
and the jet-like feature, we have verified the proper mo- 
tion of the pulsar and have shown similar motion in the 
sharp leading edge of the jet, providing strong evidence 
for its association with the pulsar. Along the major axis 
of the jet, a region of extended X-ray emission (hereafter 
counter feature), detected with ^4(t confidence, is found 
near the pulsar and opposite the jet. The jet also ap- 
pears wider, due to a significant (~5(j) region of excess 
emission, with a harder spectrum in the second epoch 
compared to the first whereas additional spectral evolu- 
tion across the major and minor axes is more ambiguous. 
Given the low confidence for many of the spectral varia- 
tions, including softening counter to the direction proper 
motion, the physical nature of the jet-like feature remains 
puzzling. Here, we briefly examine two candidate physi- 
cal scenarios for the production of the jet based on cur- 
rent results. 

To begin with, le t us first consider the model proposed 
bv lBandieral(|2008[ ). This model assumes that the jet rep- 
resents electrons that have been accelerated at the termi- 
nation shock of the pulsar wind which then leak out of the 
bow shock and diffuse into an ambient medium contain- 
ing a large-scale plane-parallel magnetic field. The ori- 
entation of the jet-like feature feature would then corre- 
spond to the direction of the magnetic field, w hich needs 
to be unusually strong (~ 45/.tG). As stated bv lBandieral 
([2008), this scenario would predict motion for the jet- 
like feature consistent with B2224-I-65, softening of the 
jet spectrum counter to the direction of proper motion 
due to synchrotron cooling of particles injected at dif- 
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ferent times and the presence of a counter feature. Our 
two epoch observations support many of these predic- 
tions though the spectral variation can not be positively 
confirmed. Additional spectral variations, such as hard- 
ening between epochs and away from the pulsar along the 
jet major axis, are difficult to explain under this scenario. 
If the jet spectrum was observed to soften away from the 
pulsar, then one may naturally expect synchrotron cool- 
ing to be responsible; however, the observed hardening 
implies acceleration within the jet feature, unlikely for 
particles confined by a magnetic field. Hardening be- 
tween epochs may suggest variation in the acceleration 
of the energetic particles, possibly due to changes in the 
pulsar wind and/or ISM properties. However, the lack of 
spectral evolution for the pulsar advocates against any 
such deviations at the terminal shock if, in fact, a sig- 
nificant portion of the pulsar emission results from the 
ISM-pulsar wind interaction. Furthermore, the relative 
faintness of the counter feature posses a challenge for the 
iBandieral ()2008[ ) model, as the energetic electrons should 
have no preference in which direction along the magnetic 
field they travel, though this issue may be rectified by in- 
ternal anisotropies and asymmetry of the pulsar wind. 

Alternatively, the B2224-I-65 extended X-ray features 
may represent the externally processed particles from 
jets, akin to those associated with AGNs. In this sce- 
nario, the main jet-like feature, counter feature and the 
Guitar Nebula represent the pulsar jets and equatorial 
wind outflow where the latter is responsible for the for- 
mation of the Guitar Nebula trailing the pulsar. X-ray 
emission is not expected in the pulsar wind forming the 
Nebula as the particles are not energetic enough to gen- 
erate detectable X-ray synchrotron radiation. The pulsar 
jets, intrinsically relativistic and containing highly ener- 
getic particles, may easily escape the bow shock but are 
ultimately confined by the large ram-pressure generated 
by the pulsar's motion. In this scenario, the spectrum of 
the jet-like feature would soften counter to the d irection 
of proper motion, similar to the IBandieral ()2008[) model, 
due to synchrotron cooling. If enhancement of the X- 
ray emission due to the high ram-pressure contributes 
significantly to the X-ray emission, then spectral varia- 
tion between epochs may arise from density variations in 
the ISM resulting in post-processing variations of the en- 
ergetic particles (i.e. through their acceleration and/or 
synchrotron cooling efficien cy). ISM variation over time 
has been seen previously by IChatteriee &: CordesI (|2004l ) 
where they found a decrease in the local density by a fac- 
tor of ^^0.7 between 1994 and 2001; however, any such 
density contrast can not be directly confirmed given the 
current two epoch Chandra observations and the rela- 
tively static pulsar spectrum. Spatial variations of the 
ISM could also explain the spectral hardening away from 
the pulsar through similar post-processing of the ener- 
getic particles. The relative faintness of the counter fea- 
ture compared to the main jet-like feature may be due to 
Doppler effects if they are indeed relativistic and aligned 
along the line of sight; specifically, if /3 cos6' ^ 0.4 where 
/3 is the bulk flow velocity and 9 is the angle along the line 
of sight. Further constraints may be possible by invoking 
more detailed jet mechanics and synchrotron theory. 

The ba llistic jet scenario was deemed unlikely by 
IBandieral (|200 8). as it would predict a ratio between 
the length and width of the jet-like feature (ijet and 



Wjet, respectively) that is too small to match observa- 
tions {Ljet ^ 0.6 pc and Wjet ^ 0.04 pc for Dkpc = 1) 
along with signiflcant bending in the jet based on com- 
parisons of the bow shock pressure, transverse momen- 
tum added to the jet through the IS M ram-pressure and 
its the initial momentum. Following iBandieral ()2008l ). we 
may estimate the expected bending angle of the jet-like 
feature under the ballistic jet model by assuming that 
the pulsar imparts a constant fraction rj of its spin-down 
energy into the jet, which has a transverse size of Wjet- 
Balancing the transverse and initial momenta, the angle 
9 describing the bending of the jet may be given by 

f^^^Q - WjetLjetsm{118°)SnmHnisMv'i: 

where 5 and are the ISM ionization fraction and 
mean molecular weight, respectively, and is the pul- 
sar proper motion. If we assume a fully ionized ISM 
composed mostly of hydrogen [5 — 1 and /i = 0.5) with 
number density njsM ~ 0.1 cm~'^ and dimensions for 
the proposed jet as taken from current observations, then 
such a ballistic jet should trail behind the pulsar regard- 
less of the value of 77; strongly suggesting against the jet 
interpretation. However, the actual jet may be very nar- 
row {Wjet ^ 0.01 pc), similar to AGN jets, and/or the 
ISM could be largely neutral in which case the expected 
bending angle would be significantly reduced. A largely 
neutral ISM would contribute little to the ram-pressure 
as the mean free path of neutral atoms is much larger 
than the width of the feature (~0.1 pc for neutral hydro- 
gen atoms with density 0.1 cm~'^); however, this would 
impose a fine-tuning problem on the jet-ISM interaction 
where the number density for the small ionized fraction 
must still be large enough to generate the confining ram- 
pressure and accelerate the energetic particles to X-ray 
emitting energies. An additional complication of the jet 
interpretation is the misalignment of the spin axis (be- 
lieved to be correlated with pulsar jets) and the direction 
of B2224-|-65's proper motion; nevertheless, the method 
by which pulsars or their progenitors acquire such high 
velocities remains unc ertain (see, for example. Hills 198g 
lYu fc Tremaindl2003^ . Gertainly, additional and more 
detailed theoretical modeling, which may require addi- 
tional observations in order to confirm apparent spec- 
tral trends, will be necessary to further test the jet and 
iBandiera (,200S) models. 

5. SUMMARY 

Using the available two epoch Chandra observations, 
we have presented an X-ray morphological and spectral 
study of the pulsar B2224-I-65 and its apparent extended 
linear X-ray-emitting features. We confirm B2224-|-65's 
proper motion over the two epochs and find little vari- 
ation in its spectrum. Additionally, we find motion in 
the sharp leading edge of the jet-like feature that is con- 
sistent with the proper motion of the pulsar. The jet 
is found to be significantly harder and wider in the sec- 
ond epoch than the first with spectral variation along its 
major and minor axis detected only marginally. These 
observations are somewhat consistent with magnetically 
confined " leakin g" of energetic particles as proposed by 
'Bandiera' (2008) as well as a typical jet interpretation; 
although there are certain limit ations to each of these 
scenarios. In the IBandieral (|2008D model, it is difficuh to 
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explain the faintncss of the counter feature as well as the 
spectral variations between epochs, if verified. On the 
other hand, it is difficult to imagine a jet being respon- 
sible for the observed emission unless the jet is actually 
very thin and unresolved, similar to AGN jet outflows, or 
if the ISM is mostly neutral, which would then impose 
fine-tuning of the ISM number density and ionization 
fraction. Detailed spectral and theoretical modeling fol- 
lowing pulsar ejecta interacting with the ISM and/or an 
ambient magnetic field, including turbulent effects within 
the ejecta, will be necessary in order to place any further 
constraints on the physical origin of the jet given cur- 



rent observations. Additional observations may then be 
required for confirmation of the observed spectral trends 
and place definite constraints on the origin of this enig- 
matic X-ray feature. 

We would like to thank the anonymous referees for 
their helpful comments on improving the manuscript and 
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